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Resolving the tension between quantum superpositions and the uniqueness of the classical 
world is a major open problem. One possibility, which is extensively explored both theoreti- 
cally and experimentally, is that quantum linearity breaks above a given scale. Theoretically, 
this possibility is predicted by collapse models. They provide quantitative information on 
where violations of the superposition principle become manifest. Here we show that the lower 
bound on the collapse parameter A, coming from the analysis of the human visual process, is 
~ 7±2 orders of magnitude stronger than the original bound, in agreement with more recent 
analysis. This implies that the collapse becomes effective with systems containing ~ 10 4 -10 5 
nucleons, and thus falls within the range of testability with present-day technology. We also 
compare the spectrum of the collapsing field with those of known cosmological fields, showing 
that a typical cosmological random field can yield an efficient wave function collapse. 



bassi@ts.infn.it 



2 



Schrodinger cat's states, namely quantum superpositions of states, each containing a large 
number of particles, are part of the jargon of quantum mechanics. Besides representing a fascinating 
challenge for experimental physics [IH3], due to their high sensitivity to external noises, they 
exemplify one of the most controversial problems in modern physics [H [5] . Quantum mechanics 
naturally allows for superpositions of macroscopic states to occur, especially in measurement-like 
processes; however, such superpositions are not part of our daily experience. 

Several ways out have been proposed. Some re-interpret or reformulate the quantum theory 
in a paradox- free way, others explicitly modify the quantum dynamics. This second possibility is 
undertaken by models of spontaneous wave function collapse [BHS] ■ 

According to these models, nonlinear and stochastic terms are added to the Schrodinger equa- 
tion. Their effect is to localize the wave function in space in accordance with the Born rule. 
Physically, the collapse is driven by a random field w(t,x) which is supposed to fill space. A 
suitably tuned coupling constant A makes sure that all quantum properties are preserved at the 
microscopic scale [9], while at the macroscopic level classicality emerges. 

As the noise w(t, x) that causes the collapse is supposed to fill space, its effects show up not 
only in laboratory experiments, but also at the cosmic scale [10J. This field interacts in fact with 
all matter in the universe. At present there is no hint whether it can be identified with any 
cosmological field considered in the literature, though some work has already been done and 
several suggested that it could have a gravitational [T2HT5] or pre-quantum [16] nature. Moreover, 
it is still not clear why it has an anti-Hermitian coupling to matter, which is necessary to ensure 
the collaps^J 

Aim of this paper is twofold: the first goal is to provide a reasonable estimate for the collapse 
parameter A. Such an estimate (its lower bound) is important, because it sets the scale at which a 
departure from the standard quantum behavior should be expected. This in turn is a crucial ingre- 
dient, in order to devise focussed experiments which would represent meaningful tests of quantum 
linearity. The second goal is to analyze the spectrum of the collapsing field as opposed to those 
of known cosmological fields. We aim at answering the question whether this field could reason- 
ably have a cosmological origin, thus whether collapse models can be considered phenomenological 
models of an underlying theory, where the collapse emerges as the consequence of the interaction 
of quantum systems with a random field of nature. 

There are several collapse models. The one that is commonly used in physical applications is the 

1 With an Hermitian coupling, one would have a standard quantum Hamiltonian with a random potential; the 
equation would be linear and no suppression of quantum superpositions would occur. 
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Continuous Spontaneous Localization (CSL) model [T7], which generalizes the original Ghirardi- 
Rimini- Weber (GRW) model [6] to quantum field theory. In the mass-proportional version of this 
model, the stochastically modified Schrodinger equation describing the dynamics is: 
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where H is the standard quantum Hamiltonian of the system, while the other terms represent 
the stochastic and nonlinear modifications which induce the collapse of the wave function. The 
operator M(x) is an averaged mass-density operator: 

M(x) = — / d 3 j/ff(x-y)Vm s aJ(y)a s (y); (2) 
m N J ^ 

in the above equation, the sum extends over particle species s of mass m s with number density 
operator at(y)a s (y), while m^v is the mass of the nucleon and g(x) is a spatial correlation function, 
conventionally chosen equal to: 

(M(x)) t = (if) t \M(x)\ipt) is the quantum average of the average mass-density operator, and Wt(x) is 
a family of standard independent Wiener processes, one for each point in space. Loosely speaking, 
the time derivative of the Wiener process: dWt(x)/dt = w(t,x), can be interpreted as a Gaussian 
stochastic field filling space, whose time correlation function is a Dirac delta. 

The fundamental (phenomenological) parameters of the model are the coupling constant 7 and 
the correlation length re, whose conventional values have been set equal to [T7j: 7 ~ 10 _30 cm 3 s~ 1 
and rc — 10 -5 cm. The collapse rate A is related to the above constant by the expression: 

The CSL model predicts that superpositions of states separated by more than the correlation 
length rc of the noise are localized with a rate T equal to [TO] : 

r = \n 2 N, (5) 

where n is the number of particles within a distance rc, and N is the number of such clusters. This 
can be seen as follows. According to Eq. 0, the decay of the off-diagonal elements of the density 
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matrix pt = E[|^)(^|] (where E[-] denotes the stochastic average) of a many-particle systems, 
which gives a measure of the collapse process, is [TT] : 

|(xV*|x") = -r(x',x")(xV t |x"), (6) 

where x' = x' 1; x^, . . . x^y (and similarly for x"). In the above equation, we have neglected the 
standard quantum evolution. The decay function T is defined as follows: 

r = I E H< - 4) + G « - 4) - 2G « - 4)] > ( 7 ) 

where the indices i,j run over the N particles of the system, and G(x) = (47rr^,) _3//2 exp[— x 2 /4r^,]. 

According to the equations above, when the particles in a superposition are displaced by a 
distance £ = |x' — x"| <C rc, their contribution to T is negligibly small. Thus, only superpositions 
with £ > rc contribute to T and trigger the collapse of the wave function. In such a case, the 
formula tells that for groups of particles separated (in each term of the superposition) by less than 
rc, the rate T increases quadratically with the number of particles while for groups of particles 
separated by more than rc it increases linearly. Thus we have the following simplified formula for 
the collapse rate: T = A n 1i where m is the number of particles of the i-th group whose mutual 
distance is less than rc- Assuming for simplicity that is constant for any i, we obtain Eq. ([5]). 
We note that the quadratic dependence of T on the number of particles — which is absent in the 
original GRW model — is a direct effect of the identity of particles. 

The conventional value of the collapse rate A corresponds to a quantum-classical threshold of 
~ 10 13 nucleons [T7]. Recently, a much stronger bound on A has been proposed [10], namely 
A ~ 2.2 x 10~ 8±2 s _1 , corresponding to a threshold of ~ 10 5 nucleons. The underlying motivation 
is that the collapse should be effective in all measurement processes, also those involving only a 
small number of particles as it happens in the process of latent image formation in photography, 
where only ~ 10 5 particles are displaced more than rc- In order for the collapse to be already 
effective at this scale, one has to increase the conventional CSL value for A by ~ 10 9=t2 orders of 
magnitude. 

Both values are compatible with known experimental data [9] • However, such a large discrepancy 
of ~ 9 orders of magnitude shows that there is no general consensus on the strength of the collapse 
process and consequently on the scale at which violations of quantum linearity can be expected to 
manifest. This does not help the setting of future experiments, testing quantum linearity. Here, 
we clarify the matter. 
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We consider the visual process, and in particular the fact that all human perceptions are clas- 
sical. Accordingly, any superposition reaching the eye must be reduced before it is transformed 
into a perception in the brain. This position has first been taken in |18j . in which authors assume 
that the reduction occurs at the very last step of the visual process, when the optical signal is 
transformed into a perception in the visual cortex. They show that by choosing A ~ 2.2 x 10 _17 s _1 
the collapse mechanism is strong enough to guarantee a classical perception. Anyway, it seems 
dubious that the collapse occurs in the latest stages of the process. Otherwise, e.g., animals with a 
simpler visual apparatus could perceive such superpositions which we consider rather unlikely. We 
thus make the assumption that the eye acts as a detector, as also proposed recently [TH]. So we 
consider the case of a superposition of few photons impinging and not impinging the retina, and 
we analyze for which particles the relative displacement in the two branches of the superposition 
is greater than rc, thus contributing to the collapse. Since we want to set a lower bound on the 
collapse parameter, we take the worst case scenario, when only ~ 6 photons, corresponding to the 
threshold of vision [20], form the input signal. When such a small number of photons triggers the 
rods in the retina, the reaction time is ~ 100 ms [21] . 

Schematically, the whole process goes as follows [22] (see Fig. 1). An incoming photon is 
absorbed by a rhodopsin in a rod disc, triggering a cis-trans transformation. In this active form, 
the rhodopsin interacts with ~ 20 transducins. The a-subunit of each transducin splits from the 
rest of the molecule, binding to a phosphodiesterase (PDE) and activating it. During this chain 
of processes these molecules diffuse on the disc surface. Next, each active PDE hydrolyzes ~ 100 
molecules of cyclic guanosine monophosphate (cGMP) to guanosine monophosphate (GMP). The 
decrease of cytoplasmic cGMP concentration causes the closure of ~ 300 ionic channels on the 
membrane of the rod [23], each of which prevents ~ 10 Na + per ms to enter the rod [24J. In this 
way the rod hyperpolarizes and a graded potential is generated: this is the electric signal, which 
is transmitted through the retina to the optic nerve. We now estimate the contribution of each 
molecule to the reduction process. 

Configurational changes of the molecules are small compared to rc and do not contribute to 
the reduction process. Only relatively large displacements count. The first relevant displacement 
is that of the a-subunits of the transducins, which remain attached to them or bind to the PDEs, 
depending on whether a rhodopsin is active or not. Since the disc surface on which the a-subunits 
can move has an area of ~ 300 r c , and reasonably assuming that the sub-units are uniformly 
distributed, then they are separated from each other by more than rc- These subunits then 
contribute with n 1 = 3.9 x 10 4 (molecular weight of each subunit) and A^ = 20, providing a factor 



6 



of n\N t = 3.0 x 10 10 to Eq. (JIJ). The next relevant displacement is that of cGMP/GMP, which 
binds to the ionic channels or diffuses in the cytoplasm, depending on its form. The concentration 
of cGMP is ~ 2^M |22| , which amounts to ~ 1 molecule in a volume of size rc, so they are 
separated by more than rc- There are ~ 2000 such molecules undergoing a hydrolysation. This 
corresponds to n 2 = 363 and N 2 = 2000, giving a contribution of n\N 2 = 2.6 x 10 8 to Eq. ([5]). 
This contribution is much smaller than the previous on^J 

The third relevant displacement is that of the Na + ions. Since the surface density of ion 
channels is 500 nmT 2 |22| . there are about 5 channels within a distance rc- We remind that ~ 300 
channels are involved in the process. In 100 ms ~ 10 3 ions enter or fail to enter the membrane 
through a channel |24j . Since each channel brings the ions inside the molecules in groups of 3 [23J , 
we distinguish two opposite situations: the more likely case case where the ions passing through a 
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FIG. 1. Visual process in the human eye. The picture shows the most relevant steps, focussing on how the 
incoming light is converted into an electric signal. 



2 Being this contribution smaller than the previous one, even considering the case that cGMP is grouped in clusters 
of size < rc, each containing ~ 100 molecules, would not affect the final result. 
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channel can be grouped in clusters of 3 molecules, distant less than rc (which means n 3 = 5 x 3 x 23, 
N 3 = 60 x 333), and the extreme case where all ions passing through a channel are distant less than 
rc (thus n 3 = 5 x 10 3 x 23, N 3 = 60). Taking into account these two possibilities, the estimated 
contribution of the third step is n^N 3 = 2.4 x 10 9 -^7.9 x 10 11 . Summing up all these contributions, 
one obtains the overall factor n 2 N, which represents the effects of one photon impinging on the 
retina. Since the threshold of vision is ~ 6 photons, we still have to multiply such a contribution 
by a factor of 6. 

When applying Eq. ^ to estimate the value of the collapse parameter, one has to decide 
when a superposition can be said to have collapsed to a localized state. We use the criterion 
adopted in |18| . according to which a superposition has collapsed when Tt ~ 10 2 , i.e. when one 
the two terms of the superposition is e 100 smaller than the other. Clearly this choice is arbitrary, 
though reasonable, and a value one order of magnitude smaller or bigger would still do the job. 
Taking this arbitrariness into account, and remembering that in our case t = 100 ms, Eq. ^ gives 
A ~ 5.0 x 10" 9±1 -h2.0 x 10~ n±1 s _1 . A further contribution comes from the fact that the potential 
travels through other cells in the retina, before entering the axon. Assuming that the number of 
particles involved in each cell is of the same order as that for the rod, and that an average of ~ 3 
cells are involved, we obtain A ~ 1.7 x 10~ 9=l:1 -j- 6.7 x 10~ 12=tl s _1 . We can summarize our result 
by saying that A ~ 10~ 10±2 s -1 . 

This value is compatible with A ~ 10 -8=t2 s _1 first proposed in |10j . Moreover, they both 
are quite surprisingly close to the prediction coming from the analysis of the Schrodinger-Newton 
equation in semiclassical gravity, according to which nonlinearities due to the interaction with the 
gravitational field should appear already for systems with ~ 10 4 nucleons [25J. At present it is not 
possible to further sharpen the numerical estimate of the lower bound. However, since A ~ 10 _8=l:2 
s _1 makes sure that, either in measurement-like processes such as latent image formation, or in the 
visual process, the collapse of the wave function becomes effective, we consider this value in the 
rest of the manuscript. If this is true, then present day technology already allows for meaningful 
tests of quantum linearity, since the collapse is predicted to be effective for systems containing 
~ 10 4 -10 5 nucleons being in a superposition^] > rc- 

A second related issue, already anticipated at the beginning of this letter, is whether the proper- 
ties of the collapsing noise are compatible with general features of cosmological fields. If nonlinear- 

3 In systems like Bose-Einstein condensates, superconductors and superfluids, a much larger number of particles 
interact coherently. However, the typical center of mass displacement is smaller than rc, thus collapse models do 
not predict a rapid loss of quantum coherence for these phenomena. Quantitative calculations show an agreement 
between the predictions of collapse models and experimental data [S]. 
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ities emerge when moving towards the macroscopic scale, and if these nonlinearities are associated 
to the interaction of physical systems with a random field filling space, then the most natural as- 
sumption is to suppose it to be a cosmological field. The first immediate question is whether such 
a field can reasonably have a cosmological nature, meaning with that, whether its most relevant 
properties (above all, its spectrum) are compatible with that of typical cosmological fields. The 
CSL model is not suitable for this purpose as the spectrum of its noise is white, comprising all 
frequencies with the same weight. Moreover, this noise represents an unrealistic thermal field with 
infinite temperature. 

Recently two generalizations of a simplified version of the CSL model (known as the QMUPL 
model |15| ) have been thoroughly investigated, the first using a finite-temperature noise [26J, the 
second using a colored noise with an exponential correlation function |27j . These models are 
particularly useful for deciding whether a cosmological fluctuating field with 'typical' properties 
can induce an efficient collapse. Equally important, they allow to follow analytically the time 
evolution of wave functions. Here we investigate how the collapse depends on the parameters of 
these models. We considered Gaussian states 



which represent typical physical states. The spread in space of the wave function, which measures 
the collapse, is: at = (2-^/of ) -1 ' 2 , the quantity af denoting the real part of at- The spread a has 
been plotted in Fig. 2 and 3. 

Fig. 2 refers to the simplified CSL model. It shows how the spread of a wave function with 
initial spread a = 5 x 10~ 7 m changes after t = 10 -2 s, as a function of the collapse strength A and 
of the number of particles. The numerical values for a and t mimic the setup of the macromolecule 
diffraction experiments [U |2S], which up to now are the interferometric experiments involving 
the largest number of particles. As the picture shows, in order for these experiments to become 
crucial tests of quantum linearity, the mass of the molecules has to be increased of ~ 2-3 orders 
of magnitude. Recent analysis [5] shows that such an improvement is already within reach. Other 
crucial future experiments involve optomechanic interferometers [3]. 

Fig. 3 (top row) shows the difference in the evolution of the spread as given by the finite- 
temperature model and by the simplified CSL model. Two temperatures have been considered: 
T = 2.73K (top row right), as for the CMBR and T = 2.73 x 1(T 3 K (top row left). In both cases the 
difference is negligible. This means that even a rather cold thermal field, according to cosmological 
standards, can induce an efficient collapse of the wave function as efficient as with the standard 
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log(cr), a =spread (to) after t = 10 2 s 




log(n), n =number of particles 

FIG. 2. Wave function spread a after a time t — 10 _2 s, as a function of particle number n and collapse 
parameter A, according to the simplified CSL model. The initial value Co — 5 x 10~ 7 m reproduces the 
geometry of the macromolecule diffraction experiments of [TJ [55] . For few particles, and for small values of 
A, the spread increases just like in the standard quantum case. Increasing A, the collapse becomes stronger. 
The black line marks where a has reduced to half its initial value, after the considered time. This can 
be taken as the threshold from quantum (no collapse) to classical (collapse) regime. In the case of the 
macromolecule diffraction experiments marked by the white strip, values of A < 10 _5 s _1 imply that collapse 
models' predictions agree with standard quantum predictions. This is consistent with the analysis of |10j 
for the CSL model. The proposed theoretical bound, represented by the white region, shows that with an 
increase of the particle number of ~ 2-3 orders of magnitude, these experiments would represent a crucial 
test for macroscopic quantum superpositions. 

CSL model. This does not mean that the two models are equivalent. On the fast collapse scale 
(1/r) they agree while on the slower relaxation scale, according to the finite-temperature model, 
the energy of any system slowly decreases and thermalizes to that of the field, while according to 
the CSL model it increases steadily. 

Fig. 3 (bottom row) shows the difference in the evolution of the spread as given by the simplified 
CSL model, and by the model with a noise having a frequency cutoff. No high-frequency cutoff 
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FIG. 3. Difference between the spread a predicted by the simplified CSL model with that given by the 
finite-temperature model (FTM), for T = 2.73K (top row left) and T = 2.73 xl(T 3 K (top row right); and 
by the colored-noise model (CNM), with cutoff at 10 10 Hz (bottom row left) and 10 2 Hz (bottom row right). 
As the color bars on the right show, the whiter the region, the greater the difference in the spreads. The 
initial spread oo and the elapsed time t are the same as in Fig. 2. At lower temperatures or lower cutoffs 
the wave function tends to collapse more slowly, which results in a bigger difference with respect to the 
simplified CSL model. Regarding the plots in the top row, the discrepancy manifests in the lower left corner 
of the plot for T = 2.73 x 10 -3 K which disappears for T = 2.73K. Regarding the plots in the bottom row, 
it is manifest in the diagonal strip for 7 = 10 2 Hz which decreases for 7 = 10 10 Hz. This diagonal feature 
exactly lies on the ridge between the quantum and the classical regime as displayed in Fig. 2. The large 
discrepancy there is due to the missing high frequencies in the noise spectrum of the colored-noise model. 
Without these high frequencies the colored-noise model is not able to reproduce the sharpness of the ridge 
predicted by the simplified CSL model. 



affects the collapsing properties of the model in an appreciable way. These results can be compared 
with the behavior of typical cosmological fields such as the CMBR, the relic neutrino background 
and the relic gravitational background. The spectrum of the first two have a cutoff (measured or 
expected) at ~ 10 11 Hz, while the spectrum of the third one probably lies at ~ 10 10 Hz [29J. All 
these cutoffs as well as that of ~ 10 15 Hz proposed in [8] for w(t, x) ensure a rapid collapse of 
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the wave function. While the collapse is robust over a large range of cutoffs, other effects, such 
as the emission of radiation from charged particles, highly depend on the spectrum of the noise 
correlator [5U] . 

The message which can be drawn is that a cosmological field with 'typical' properties can induce 
an efficient collapse. A great challenge is to test the existence of such a field. Our analysis, like 
that of [10], points towards a much stronger collapse effect than previously thought. This means 
that crucial experiments testing quantum linearity are within reach. 
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Appendix 1: Finite temperature model. According to this model, the time evolution of at 
is given by |26j : 
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a + ib tanh ( bt + k 



(9) 



where k sets the initial condition a>o, while 
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a\ 2 + 2i — -n 2 . (10) 



2k B T V h 

In the above expressions, T is the temperature of the random field, mo is the nucleon's mass and 
k,B is Boltzmann constant. By taking the infinite temperature limit (i.e. a — > 0), one recovers the 
time evolution as given by simplified CSL model [15] . 

In using the formulas above we have replaced n with n 2 . This way we take into account that the 
collapse rate of the CSL model depends quadratically on the number of particles whenever they 
are bounded in systems of size smaller than rc- This in particular is the case for macromolecules. 

Appendix 2: Colored noise model. Assuming an exponential correlation function D(t,s) for 
the noise [27J: 

D(t,s) = ( 7 /2)e-^- s \, (11) 
where 7 gives the cutoff, the function at can be written as: at = At — (£>f /4(ao +«4t))j where we 
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have defined 



+ (dk sinh t^i — c cosh t^i) sinh t^i + , (12) 



Bt = ~ "^0~ Z_> Vk 1°* cosh U S* + b ~k smh u k* 

+ dj, cosh ffct — c sinh t^t] , (13) 
6( = ^ [2c + (aj. coshfgt + 6^ sinhu^t) smli-u^i 

k 

+ (dk sinh vj.t — c cosh t^i) cosh v^t] (14) 
with /c £ {+, — } and A; denotes the opposite sign. The new coefficients are a± = 7t;±[i;± =F (}, 



b± = f±[i>i T 7 2 C]> c = v^_v\, d± = — 7t> with: 



Also in this case, in using the above formulas we have replaced n with n 2 to reproduce the fact 
that the collapse rate of the CSL model depends quadratically on the number of particles. 
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